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Introduction 


This report summarizes the work done in the period January 1994 through June 1994 on 
the "Modeling of Outgassing and Matrix Decomposition in Carbon-Phenolic Composites" 
program. 


Progress 

Two threads of research have been followed. First, the thermodynamics approach was 
used to model the chemical and mechanical responses of composites exposed to high temperatures. 
This approach was first proposed by Sullivan [1], The thermodynamics approach lends itself 
easily to the usage of variational principles. This thermodynamic-variational approach has been 
applied to the transpiration cooling problem. In transpiration cooling, a porous solid is cooled by 
gas flowing from the cool end to the hot end. The governing equation for the temperature 
distribution was obtained for the steady- state case. 

The second thread is the development of a better algorithm to solve the governing equations 
resulted from modeling. The resulting nonlinear, partial differential equations are vety complex. 
The complexities are increased even more when detailed material models are included such as 
multi-stage chemical reactions and non-ideal gases. Previously, implicit finite difference method 
was used to solve the governing equations. Although the implicit method is relatively stable, it 
was found that detailed material models cannot be easily incorporated into the implicit finite 
difference method. Explicit finite difference method is explored as an alternate numerical method 
to solve the governing equations. In the explicit method, detailed material models can be included 
into the algorithm much easier than the implicit method. The explicit method is also highly suitable 
for solution on massively parallel supercomputers such as SCOUT, MITs new ARPA-sponsored 
Thinking Machines CM-5. 


To demonstrate the feasibility of the explicit scheme in solving nonlinear partial differential 
equations, a sample problem was solved using SCOUT. The sample problem was a transpiration 
cooling problem, which resembles the composite ablative problem in that the solution requires 
solving simultaneously the continuity, energy, and momentum equations. Some interesting 
transient behaviors were captured such as stress waves and small spatial oscillations of transient 
pressure distribution. The results of this study were presented at NASA-Langley in April 1994. A 
set of viewgraphs from this presentation are attached. 




Current Status 


The thermodynamic-variational approach is currendy been pursued. Using the approach, 
governing equations for material responses of composites exposed to high temperatures will be 
derived. Biot has done some extensive work in developing the thermodynamic-variational 
approach in references [2,3]. Our work will follow Biot's development closely. 

The ablative composite code CHAR is being re-coded using the explicit finite difference 
method on the CM5 Massively Parallel Machine that MIT has acquired recently. After CHAR- 
CM5 is completed, the new thermodynamically-based material and reaction models will be 
included into the code to perform parametric studies. 
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COMPOSITE ABLATORS 



PI .V LIFTING 

Very large heat loads require ablative materials 


Modeling Material Response 
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Modeling - Approach 
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Quantitatively Accurate Model Allows: 
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Reaction Zones 
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Pyrolysis Zone and Moisture Evaporation Zone are very small 
Large pressure developed just below the Pyrolysis Zone 
Ply-lift is observed to occur just below the Pyrolysis Zone 


Internal Pressure 






Composite Plate 




Insulated, Impermeable Base 


Control 

Volume 



Conservation of Mass 
Conservation of Energy 
Equations of State 
Darcy's Flow Law 
Reaction Kinetics 



THERMOCHEMICAL MODEL 











Immediately Before Damage 
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Stress Ratio, R 




Moisture Loss 
MC 0 - MC (percent) 


Pressure Difference 
*P = P-P b (MPa, 
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Current Research 

Qualitative predictive capability achieved 
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ADVANCED NUMERICAL METHODS 
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- More difficult to program 

- Nonlinear equations need large iterations for convergence 




Unknown at t+At expressed in terms of knowns 

- Simple algorithm 

- Need small time steps 



EXAMPLE PROBLEM: 
TRANSPIRATION COOLING 
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Check Case #2 
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Parametric Studies 
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Variation in Pi 
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Variation in Permeability 
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V ariation in Permeability 
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Variation in Porosity 
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Variation in Porosity 



Thickness (cm) 




Examples 




a c n 



CA 




eo m 
cs co 


ii ii » » 

cCi£h£ 


permeability = 100 x 10 1S m 2 
viscosity = 1 x 10 5 kg/m sec 

porosity = 0.05 



Transient Temperature 
Distributions 
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FUNCTIONALLY GRADED 

MATERIALS 




ALTERNATE MATERIAL DESIGNS 


MULTIPLE MATERIALS 


Design Parameters 

Material Choice 
(Two or more) 


Constraints Problem 


Available Materials 
Compatibility 
Joining Techniques 


Stress Concentrations 
Weight 
Expense 




COMPOSITES 


Design Parameters 

Constituents 
Volume Fractions 
Geometry 


Constraints 

Available Materials 
Compatibility 
Processing Techniques 


• Problems 

Spatial Uniformity 
3rd Dimension 
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FGM MATERIAL DESIGNS 


Hot Environment 



Design Parameters Constraints Problem s 

Constituents Available Materials TBD 

Volume Fractions Compatibility 

Geometry PROCESSING 

ALL FUNCTIONS 
OF POSITION 




DESIGNING WITH FGM 

♦ FGM frees up the spatial distribution of material properties as design 
parameters 6 


• This is a key advantage in problems with multiple functional 
requirements 


EXAMPLE 


Hot Environment 




Mechanical 

Loading 


Cool Interior 


TRADITIONAL MATERIAL CHOICE 


Desien Parameters 
Material Choice 


Constraints 
Available Materials 


Problems 

All props set by single 
choice 

May be no solution 
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INTEGRATED FGM DESIGN 


TRADITIONAL MATERIAL SELECTION PROCESS 


NEW MATERIAL 


STRUCTURAL 
DESIGN AND 
MATERIAL 

CUT PPTrnKT 

DEVELOPMENT 

PICK 
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NEEDS DRIVE NEW DEVELOPMENT 
TIME CONSTANT = 10 YEARS 


Decoupled material and structural design 

Very long time before new materials are ready for application 




INTEGRATED FGM DESIGN 


Fully integrated design, and dose cooperation between designers 
analysts, and processors, critical for rapid exploitation of FGM 


NEEDS 

DRIVE 

NEW 

DEVELOPMENT 
T = YEARS 



CLOSE 

COOPERATION 
T = MONTHS 


Parallels major initiative in the MIT Aeronautics and Astronautics 
Department in Engineered Materials 

Integration of material design into aerospace structural engine ring 
• Research and educational initiatives ° 






TOOLS FOR INTEGRATED DESIGN 


PROCESS MODELING 


Given process parameters, what is the microstructure as a function of 
position? 


MICROMECHANICS 

Given constituent properties, volume fractions, architectures and 
microstructures, what are the local bulk properties? 

GRADIENT STRUCTURE ANALYSIS 


Given local bulk properties as functions of position, what ar the 
distributions of temperature, stresses, etc., also as functions f 
position? 

GLOBAL STRUCTURAL ANALYSIS 


What is the response of the entire structure? 


1 793’C 




All these analyses must be coupled 



EXAMPLE DESIGN ANALYSIS 

Computer code adapted from charring ablator code CHAR 
Coupled thermal, porous flow, chemical reaction, and thermo-poro- 
elastic structural solutions 

Transient solutions for arbitrary thermal BCs, and several structural 
cases (coatings, plates) 

Simple FGM input 
Crude micromechanics 

Limited material failure and damage modeling capabilities 


Heat Flux + Pressure 



Cooling gas 

User-Specified Graded Material 



EXAMPLE CASE 


Environment = 500°C 



Cooling gas = 0°C 
No Transpiration 


Silica/Copper with Various Gradients 
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EXAMPLE CASE RESULTS 




0.03 





SUMMARY AND CHARGE 

To apply FGM concepts to practical design problems as quickly as possible 

• Bottleneck disciplines must be identified and worked on 

• Existing processing, analysis and design capabilities must be i legated 
for FGM application 

It is our hope that this workshop will 

• Identify what can and cannot be done with existing technique 

• Aid in the integration of existing knowledge 

• Identify critical areas for research and development 
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several proressors will supervise each student 



HEAT TRANSFER PROBLEM 
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